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SUMMARY 


An  approximate  method  for  determining  the  far  field  flow  properties  of  a 
rocket  exhausting  into  a  vacuum  is  presented.  The  technique  described  is  based  on 
radial  flow  and  assumes  a  density  variation  both  along  and  normal  to  the  jet  center- 
line.  The  resulting  jet  flow  field  is  compared  with  solutions  by  the  method  of 
characteristics  and  is  shown  to  have  good  agreement.  Example  calculations  are 
provided  to  demonstrate  the  present  method. 
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NOMENCLATURE 


a  Speed  of  sound 

A  Area 

E  Constant  defined  by  Equation  (13) 

E'  Constant  defined  by  Equation  (26) 

h  Radius  from  nozzle  exit 

M  Local  Mach  number  at  the  point  (h,  6  ) 

M'  Local  Mach  number  at  the  point  (h,  Q  )  in  the  flow  field  of  a  sonic  nozzle 

N  Constant  defined  in  Equation  (4) 

P  Pressure 

r  Radius 

S  Control  surfa  e 

V  Velocity 

X  Axial  distance  from  nozzle  exit 

Y  Distance  normal  to  jet  centerline 

7  Gas  ratio  of  specific  heats 

p  Gas  density 

Subscripts 

e  Nozzle  exit  conditions 

*  Nozzle  throat  conditions 

o  Stagnation  (chamber)  conditions 

s  Refers  to  spherical  cap 
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OBJECTIVES 


The  objective  of  this  technical  note  is  to  describe  an  approximate  method 
for  determining  the  far  field  flow  properties  for  a  rocket  exhausting  into  a  vacuum. 


CONCLUSIONS  AND  RECOMMENDATIONS 


An  approximate  method  for  determining  the  far  field  flow  properties  for  a 
rocket  exhausting  into  a  vacuum  has  been  presented.  A  comparison  of  Mach  number 
distribution  both  along  and  normal  to  the  jet  axis  with  solutions  by  the  method  of 
characteristics  indicates  good  agreement.  The  present  method  should  provide  an 
adequate  approximation  for  rocket  flow  field  problems  in  a  vacuum.  Recommendations 
are  the  following: 

(1)  Investigate  further  the  limits  of  this  approximate  solution  for: 

(a)  Nonvacuum  conditions. 

(b)  Nozzles  with  noncircular  cross  sections  at  the  exit  plane. 

(2)  Develop  an  Improved  approximate  method  which  describes  the  exhaust 
flow  past  90  degrees.  (The  present  method  is  limited  to  expansion 
angles  less  than  90  degrees.) 
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SECTION  I 


INTRODUCTION 


When  a  rocket  exhausts  into  a  vacuum,  the  flow  undergoes  a  large  expan¬ 
sion  and  the  resulting  free  jet  of  gas  occupies  a  considerable  volume  in  space. 

The  presence  of  this  large  free  jet  of  gas  can  produce  serious  aerodynamic,  aero- 
thermal,  and  communications  problems.  Some  of  these  problems  are:  the  effects  of 
Jet  impingement  forces  on  the  dynamics  of  nearby  surfaces;  the  exhaust  blast  on 
the  surface  of  the  moon  during  ascent  or  descent  of  a  rocket;  the  heating  of  nearby 
surfaces  submerged  in  the  rocket  plume;  and  the  effects  of  the  exhaust  gas  on 
visual  or  other  guidance  information. 

In  order  to  investigate  these  problems,  an  adequate  definition  of  the 
rocket  exhaust  flow  field  is  necessary.  In  most  ci.ses,  the  method  of  character¬ 
istics  has  been  utilized  to  determine  properties  in  the  rocket  flow  field.  How¬ 
ever,  the  procedure  is  quite  laborious,  costly,  and  generally  requires  machine 
computations.  To  provide  a  more  convenient  and  practicable  means  for  defining  the 
rocket  flow  field  in  a  vacuum,  an  approximate  method  is  described  in  this  technical 
note.  The  technique  can  be  utilized  to  calculate  pressures,  densities,  temperatures 
and  Mach  numbers  in  the  rocket  flow  field  from  a  knowledge  of  only  the  nozzle  geom¬ 
etry  and  the  ratio  of  specific  heats  for  the  gas. 
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SECTION  2 


ANALYSIS 


Consider  the  exhaust  Jet  of  an  ideal  inviscid  gas  issuing  isentroplcal  1 y 
from  a  nearly  parallel  nozzle  (of  exit  radius  rg  and  exit  Mach  number  Mg)  Into  a 
vacuum.  The  flow  Is  continuous  and  steady,  and  the  ratio  of  specific  heats,  7  » 
is  ass imed  constant. 


At  distances  large  compared  to  the  nozzle  dimensions,  theoretical  solu- 
tions  for  the  exhaust  flow  in  a  vacuum  show  that  the  flow  field  approaches  radial 
flow,  i.e.,  the  streamlines  are  straight  and  appear  to  emanate  from  a  common  nource 
(see  Figure  1).  The  basic  solution  of  radial  flow  stipulates  that  the  mass  flux, 
pV,  varies  as  1/h^,  where  h  is  the  radial  distance  from  the  nozzle.  In  a  vacuum, 
the  pressure  rapidly  approaches  zero  so  that  the  velocity  can  be  assumed  constant. 
In  this  case,  the  density  takes  the  form  /  ~  1  / h ^ . 

Most  of  the  mass  and  momentum  of  the  jet  are  concentrated  near  the  jet 
centerline,  so  that  the  variation  of  density  in  the  direction  normal  to  the  jet 
centerline  must  be  considered.  Reference  1  assumes  that  the  density  on  a  spherical 
cop  at  a  distance  h  from  the  nozzle  exit  is  of  the  form: 


(cos  0 ) 


K 


(1) 


where  '■>  is  tne  azimuthal  angle  measured  from  the  jet  centerline,  and  B  and  K  are 
constants.  By  integrating  the  conservation  of  mass  and  momentum  equations,  and 
neglecting  terms  of  order  1/Mg,  Reference  1  has  obtained  a  first-order  solution  for 
the  density  profile: 


JL 

P„ 


K 

2 


.  cos 


or 


(2) 
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where  K  * 


e 


1  (  7  -1)  M* 


The  above  expression  has  t 


restriction  that  M  » 

e 


and 


In  the  present  analysis,  a  more  general  solution  which  includes  the  effects 
of  higher  order  terms  is  sought,  l.e.,  it  is  desirable  to  remove  the  restriction  in 
Equation  (2)  that  Mg  »  1.  The  basic  approach  to  the  problem  is  that  of  Reference  1. 
A  density  profile  similar  in  form  to  Equation  (1)  is  assumed.  The  continuity  and 
momentum  equations  are  then  solved  simultaneously  to  obtain  the  density  distribution 
in  the  jet  flow  field.  From  this  result,  a  generalized  method  for  calculating  con¬ 
tours  of  constant  Mach  number  for  a  rocket  flow  field  in  a  vacuum  is  derived. 


The  system  under  analysis  is  shown  in  Figure  2  (a)  and  (b) .  The  flow  is 
considered  radial  and  axially  symmetric.  The  surface  S  is  the  sum  of  the  hemispher¬ 
ical  surface  (Ag)  of  radius  h,  and  its  projection  on  the  Y-Z  plane,  Ay_£.  Flow 
enters  the  control  surface  only  through  the  nozzle  exit  area,  Ag,  and  exits  through 
the  surface,  Ag. 


is,* 


On  the  spherical  cap,  the  velocity  is  assumed  radial  and  constant.  That 


V 


8 


V  tl 
max 


(3) 


where  Vmax  is  the  maximum  velocity  obtainable  by  expanding 
ia  a  unit  vector*  everywhere  normal  to  the  surface  S.  The 
surface,  Ag,  is  assumed  to  be  of  the  form: 


to  zero  pressure,  and  n 
density  profile  on  the 


(cos  6 


(A) 


where  E  and  N  are  constants. 


For  steady  flow,  the  integral  form  of  the  continuity  equation  is: 


(V  •  n)  dS 


0 


Over  the  surface  S  (see  Figure  2),  the  above  integral  is  written: 


A, 


p  (V  •  r)  dA  -  (  p  V  A  ) 

!  s  8  r  e  e  e 


(5) 


Introducing  the  expressions  for  V  from  Equation  (3)  and  p  from  Equation  (A),  then 

s 


♦Vector  quantities  are  symbolized  by  bold-face  letters. 
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Equation  (5)  becomes: 

7T/2 


1  f  1 

1  v  ) 

1 

1 

e 

h_ 

V 

r 

1  ] 

i  max| 

i  e 

(cos 6  ) 

Integrating  and  solving  for  E,  we  obtain: 

E  =  N+l 


N 


V  n  •  n 

max 


2  7T  h  si n© 


(6) 


J  0  «  l  V  A 

re  e  e 


(conti nui ty) 


(7) 


Conservation  of  momentum  need  be  considered  only  in  the  X-direction,  since 
the  flow  is  axially  symmetric.  The  integral  form  of  the  conservation  of  momentum 
equation  in  the  X-direction  is: 


ffpa 


i p(V  •  n)  dS  (V  •  i)  +  £j)p(*i  •  I)  dS  -  0 

S'  S 

For  the  surface  S,  the  above  integral  becomes: 


j  f  |"p(Vg  •  n)  (Vs  •  I)  +  p(ti  •  1) 


dA  =  p  V  +  p  I  A 
8  |  ee  Ke|  e 


(8) 


Substitution  of  the  expressions  for  Vg  and  p  (from  Equations  (3)  and  (4) 
respectively)  into  Equation  (8)  results  in  the  integral: 

7T  /  2 


V  /  '2 

,  E  I  e  h  ,  .  ,N 

/e  2  ”  17“  (cofl  0) 

max  e 


!  V  ti  •  ti 
max 


V  ti  •  t 
max 


2  7T  h  sin© 


-  p  V*  +  p 

'  o  o  c 


e  e 


■  P  V‘ 
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1  +7iil 
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d0 


Solving  Equation  (9)  for  E,  we  obtain: 

E  *  (N  +  2) 


V  \ 
e 


V 

^  max  j 


1  +  7  M2 


(momentum) 


(10) 


The  constant  E  can  now  be  determined  by  solving  Equations  (7)  and  (10) 
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simultaneously.  Thus, 


(ID 


The  velocity  ratio,  (V  /V  ),  is  related  to  the  nozzle  exit  Mach  number,  M  ,  by 
the  adiabatic  relationship3* 


max  1,1+ 


(7  -dm; 


(12) 


Then,  Equation  (11)  becomes: 


1  + 


7m; 


E  =  — 


v  + 


field  is  obtained 


Jl 
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<7 -DM* 

Equation 

(4), 
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l  e 

fn_\ 
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'»  +  “T 

\  ^Me 


(13) 
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(E-l) 


cos 


0) 


(14) 


where  E  is  defined  by  Equation  (13)  and  ( Vg /V^^)  by  Equation  (12),  It  is  seen  that 


both  E  and  (V  /V  )  are  functions  of  only  M  and  the  ratio  of  specific  heats,  7  . 
e  max  e 

Thus,  for  a  given  gas  and  nozzle  exit  Mach  number,  M  ,  the  density  depends  only  on 
the  location  in  the  flow  field  (h,fr). 


Assuming  that  the  flow  inside  of  the  nozzle  is  isentropic,  the  Mach 
number  distribution  in  the  flow  field  can  be  derived  from  thr  isentropic  relation¬ 
ship  between  density  and  Mach  number. 
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That  is, 


£_ 

Pn 


=  ^1  +  M2 


1-7 


where  p  is  the  stagnation  (chamber)  density. 


Solving  the  above  expression  for  the  Mach  number,  we  obtain 

^  1/2 
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<  >  ! 

LeJ 
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|7  -1 

p°l 

The  density  ratio,  ,  can  be  expressed 


(15) 


(16) 
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Substituting  Equation  (14)  into  Equation  (17),  then 


e  ,\V h  ,-2 


P. 


2  V  '  ,  r 
'  max  ,  i  e ; 
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(E-l) 
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From  continuity  considerations, 


P  v 
e  e 


Thus,  Equation  (18)  becomes: 
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Introducing  the  isentropic  relationships. 
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P* 

Po 


~y  +  1 
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+1\  l'7 


(17) 


(18) 


(19) 


(20) 


(21) 
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(22) 


max 


1  1/2 


IJLdj 

,v+1  I 


Then  Equation  (20)  is  written: 


,  -2 
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(23) 


Combining  Equation  (23)  with  Equation  (16),  the  final  form  for  the  Mach  number  is 
obtained: 


M  = 


1 


2  J  '  E 
7-1  \  2  !  A, 


f2 

j  (cos  9 ) 


(E-l) 
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(24) 


1/2 


where  E  is  given  by  Equation  (13). 


Significant  to  note,  for  a  given  gas  and  nozzle  geometry,  the  Mach  number 
depends  only  on  the  location  in  the  flow  field  (h,  0).  The  above  expression, 
Equation  (24)  is  general  and  can  be  used  to  find  local  Mach  numbers  in  the  flow 
field  for  any  specified  M  and  7  .  The  main  restriction,  implied  in  the  radial  flow 

assumption,  is  that  —  »1 . 


Equation  (24)  has  been  applied  to  construct  contours  of  constant  Mach 
numbers  in  Figure  3  for  Me  =  2.21  and  7  =  1.24.  Axial  distances  from  the  nozzle 
exit  are  given  by  X  and  distances  normal  to  the  jet  centerline  by  Y.  Evidently 
h^  *  X^  +  Y^,  and  0  =  Tan~l  Y/X.  The  distances  X  and  Y  are  normalized  by  the 
nozzle  exit  radius,  re.  A  numerical  example  of  Equation  (24)  is  also  provided  in 
the  Appendix,  Example  1. 


Examination  of  Equation  (24)  suggests  that  for  each  set  of  M(  and  7  ,  a 
separate  Mach  number  chart  like  Figure  3  must  be  constructed.  For  many  variations 
in  Me  and  7  ,  the  computations  become  quite  excessive.  However,  by  considering  a 
specified  gas  (7  fixed),  a  mathematical  simplification  of  Equation  (24)  is 
possible.  Let  M’  designate  die  local  Mach  number  at  the  point  (h, 0 )  in  the  flow 
field  of  a  sonic  nozzle  (Me  *  1.0);  and,  M,  the  local  Mach  number  at  the  same  point 
(h,  0 )  in  the  flow  field  of  a  nozzle  with  an  exit  Mach  number,  Mg,  other  than  one. 
The  ratio  of  specific  heats,  7  ,  is  the  same  for  both  cases. 


In  the  general  case,  for  a  nozzle  exit  Mach  number  other  than  one,  the 
local  Mach  number,  M,  is  given  directly  by  Equation  (24).  Combining  Equations  (25) 
and  (24), 


(27) 


For  —  »1,  the  terms  of  the  left  hand  side  of  the  numerator  of  Equation  (27)  are 
e 

large  compared  to  unity.  That  is, 
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By  a  similar  argument,  the  denominator  of  Equation  (27)  can  be  written: 
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From  Equations  (28)  and  (29),  Equation  (27)  now  becomes: 
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(  cos  b  ) 


l-  y 
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(30) 


where  E  is  defined  by  Equation  (13)  and  E'  by  Equation  (26). 


Having  once  constructed  the  flow  field  for  a  sonic  nozzle  from  Equation 
(25)  for  a  given  gas,  the  Mach  number,  M,  in  the  flow  field  of  a  nozzle  with  exit 
Mach  numbers  other  than  one  can  be  readily  determined  from  Equation  (30).  In  other 
words,  for  a  given  gas  (7  fixed),  plots  of  Equation  (25)  and  (30)  are  sufficient  to 
establish  the  flow  field  Mach  numbers  for  several  nozzle  exit  conditions. 

Equations  (25)  and  (30)  are  plotted  in  Figures  9(a)  and  4(b),  respectively,  for 
7  *1.4.  An  example  calculation  to  demonstrate  the  use  of  these  two  charts  is 

provided  in  the  Appendix,  Example  2. 


In  preliminary  des'gn  and  proposal  studies,  wherein  various  nozzle  exit 
conditions  and  gases  may  be  ur  er  investigation,  plots  of  Equations  (25)  and  (30) 
similar  to  that  in  Figures  4(a)  and  4(b)  for  other  gases  should  be  of  considerable 
value . 
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SECTION  3 


RESULTS  AND  DISCUSSION 


Flow  field  Mach  numbers  computed  by  the  present  approximate  method  are 
compared  with  available  solutions  by  the  method  of  characteristics  for  sever/  1  gases 
and  nozzles  exhausting  into  a  vacuum. 

Comparisons  are  presented  in  Figures  5  and  6  for  h/r(  <  35.  Figure  5 
shows  a  comparison  over  the  entire  flow  field  for  Mg  *  2.94  and  7*1-4 
(Reference  2,  for  a  15 -degree  conical  nozzle);*  and,  in  Figure  6,  the  Mach  number 
distribution  along  the  jet  centerline  is  compared  for  Mg  *  1.0  and  7  *1.4 
(Reference  3,  for  a  nearly  sonic  orifice) .  Also  shown  in  Figure  6  is  the  first 
order  solution  based  on  Equation  (2) . 

Figures  7  and  8  show  Mach  number  comparisons  for  higher  values  of  h/ rg 
(between  40  and  600).  The  entire  flow  field  is  compared  in  Figure  7  for  Mg  *  3.93 

and  7  =  1.26  (Reference  4  for  a  15-degree  semidivergence  nozzle  angle).  Axial  Mach 
number  comparisons  are  presented  in  Figure  8  for  Me  =  4.0  and  7  =  1.24  (Reference  4 
for  a  20-degree  semidivergence  nozzle  angle)  . 

Examining  the  axial  Mach  number  distribution  in  Figures  5  through  8  (shown 
directly  in  Figures  6  and  8  aid  along  the  axis  0  *  0  in  Figures  5  and  7)  the 
present  method,  in  general,  is  seen  to  overestimate  the  Mach  number  (relative  to  the 
Mach  number  by  the  method  of  characteristics)  at  a  given  value  of  X/re  along  the  axis. 
Inasmuch  as  the  comparisons  presented  in  Figures  5  and  6  are  for  relatively  low 
values  of  h/re,  the  agreement  with  characteristics  solutions  is  surprisingly  good. 
Note,  for  example,  in  Figure  5  at  X/rg  =  25.4,  and  in  Figure  6  at  X/ rg  ■  20,  the 
axial  Mach  number  computed  by  the  present  method  are  in  error  by  roughly  2  percent, 


*  In  the  development  of  this  approximate  method  the  mass  flux  /  V  A  into  the 
control  surface  was  taken  to  be  constant,  i.e.,  parallel  nozzle6  6For  a  conical 
nozzle  this  is  still  a  reasonable  approximation  and  the  method  should  be  valid 
for  conical  nozzles  as  well. 
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relative  to  these  by  the  method  of  characteristics.  At  larger  values  of  h/rg,  where 
the  radial  flow  assumption  is  more  appropriate,  better  agreement  is  evidenced.  In 
Figure  7,  at  X/'ie  *  134,  and  in  Figure  8  at  X/re  *  500,  the  relative  errors  in  axial 
Mach  number  are  on  the  order  of  1  percent.  An  examination  of  the  axial  Mach  number 
distribution  in  Figure  6  reveals  that  the  present  method  represents  a  substantial 
improvement  over  the  first  order  solution  (Reference  1).  At  X/rg  *  20,  the  relative 
error  in  axial  Mach  number  is  reduced  from  16  percent  (first  order  solution)  to 
2  percent  by  the  present  solution.  Note  also  that  the  present  solution  converges 
toward  the  characteristics  solution  as  X/rg  increases;  whereas,  the  first  order 
solution  does  not  show  this  trend. 

Normal  to  the  jet  centerline,  the  present  method  predicts  the  general 
shape  of  the  Mach  contours  fairly  accurately  (see  Figures  5  and  7).  For  the  same 
constant  Mach  number,  it  is  seen  that  the  two  curves  (approximate  and  characteris¬ 
tics  solutions)  cross  each  other.  Evidently,  at  the  point  of  intersection,  the 
agreement  between  the  two  solutions  is  exact.  Over  the  entire  flow  field,  however, 
no  consistency  is  evidenced  with  respect  to  the  angle  9  at  which  this  crossover 
occurs  . 


The  deviations  discussed  above  are  attributed  to  the  assumed  density  pro¬ 
file  on  which  the  present  method  is  based.  The  assumed  density  profile,  Equation  (-'•)» 
does  not  satisfy  the  boundary  conditions  at  the  edge  of  the  jet.  That  is,  the 
cosine  distribution  of  Equation  (4)  presumes  a  boundary  (M  =  00  )  at  9  *  90  degrees, 
whereas,  in  theory,  the  exhaust  flow  from  a  nozzle  into  a  vacuum  can  expand  to  an 
angle  greater  than  90  degrees,  depending  on  the  gas  ratio  of  specific  heats,  and  the 
nozzle  geometry.  Better  correlation  with  the  method  of  characteristics  should  be 
realized  by  considering  the  flow  past  90  degrees.  Evidently,  a  oensity  distribution 
other  than  the  present  cosine  function  must  be  assumed. 

No  general  criteria,  applicable  for  all  gases  and  nozzle  geometries,  has 
been  established  with  regard  to  the  lower  limit  (^  )  of  applicability  for  the 

*e 

present  method.  However,  an  examination  of  several  theoretical  streamline  patterns 
for  nozzles  exhausting  into  a  vacuum  (References  2  and  5)  has  indicated  that  for 
7  *1.4  and  Me  3,  the  lower  limit  may  be  taken  as  h/ re  =  10 

The  present  approximate  method  can  be  extended  to  estimate  Mach  number 
contours  for  a  rocket  exhausting  into  a  near  vacuum  environment  (such  as  in  a 
vacuum  chamber).  Whereas  the  flew  field  in  a  vacuum  is  shock  free,  the  flow  behind 
a  rocket  exhausting  into  a  finite  pressure  region  (at  rest)  is  characterized  by  a 
curved  shock  which  closely  follows  the  jet  boundary  and  eventually  returns  to  the 
Jet  axis  as  a  normal  shock  further  downstream  of  the  nozzle  exit.  The  adjustment  of 
flow  to  the  ambient  pressure  occurs  in  the  outer  region  between  the  shock  and  jet 
boundary.  Interior  of  the  shock,  the  flow  is  unaffected  by  the  ambient  pressure  and 
is  identical  to  the  corresponding  portion  of  a  vacuum  jet  flow  field  having  the  same 
Me  and  7  • 
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APPENDIX 


EXAMPLE  CALCULATIONS 


The  present  method  can  be  utilized  to  determine  properties  at  the  point  P 
in  the  flow  of  a  nozzle  exhausting  into  a  vacuum.  Tne  method  is  applicable  for 
h/re  »1.  For  a  given  gas  specific  heat  ratio  (7),  nozzle  dimension  (re),  and 
either  (Ae/A*)  or  Me,  properties  at  the  point  P  can  be  determined  by  either  of  the 
following  two  methods: 

(1)  The  basic  solution  derived  from  continuity  and  momentum 
considerations:  Equation  (24)  can  be  solved  directly  for 
the  Mach  number. 

(2)  The  generalized  solution  based  on  a  sonic  nozzle  flow 
field.  Combined  plots  of  Equations  (25)  and  (30)  provide 
a  rapid  means  for  calculating  the  flow  field  behind  a 
nozzle  having  exit  Mach  numbers  other  than  one. 

Example  calculations  are  presented  for  7  -1.4.  In  Example  1,  the  local  Mach 
number  at  the  point  P  is  determined  directly  from  Equation  (24).  The  same  problem 
is  solved  in  Example  2  using  the  generalized  procedure  based  on  a  sonic  nozzle  flow 
field  (Equations  25  and  30) . 
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EXAMPLE  CALCULATIONS 


Example  I  (Direct  application  of  Equation  24) 

Given:  Me  =  3.0,  rg  =  2  inch,  7=1-4 

Calculate:  Local  Mach  number,  M,  at  the  point  h 

h  A* 

(a)  Determine  —  ,  E,  and  — —  : 

r  „  A 

e  e 


40  inch,  G 


— rrr  from  NACA  1135  for  M  -  3.0,  7-1.4 

4.235  e 


40  in . 
2  in. 


<7 -dm"  i  7m; 


Equation  (13) 


(1 -4) (3) 

V  +  - ~ 2 

V  (1 .4-1) (3) 
6.43 


1  .4(3  ) 


(b)  The  Mach  number,  M,  is  calculated  from  Equation  24: 


_j_7  e  h_\  . 

2  ~  r  1  <cos0> 

l-  ‘I  ei 


‘•‘ixt I  l-y  1^ II  1/2 
[  2  j  7  +1 


2  6.43 

1.4-1  )  2 


(20)  2  (cos  30°) 


o/6  43-1)|l^±I\  ^ 

I  2  I 


'l.4-l\  1/2 
1.4+lj 


(1-1.4) 


1  i  .78 
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EXAMPLE  CALCULATIONS 


Example  2  (From  combined  plots  of  Equations  25  and  30) 

Given:  Me  ■  3.0,  re  *  2  inch,  7  “  1 .4 

Calculate:  Local  Mach  number,  M,  at  the  point  h  =  40  inch  and  9  *  30° 


*  h 

(a)  Determine  — ■—  and  — 
A  r 

e  * 


- — r—  from  NACA  1135  for  M  =  3.0,  7=14 
4  .235  e 


h_ 
r , 


h_  ->  l _ e 

V  A, 


V  4.2  35 

2 


h 

r. 


41  .2 


(b)  Find  M'  from  Figure  3(a)  (sonic  nozzle  flow  field,  7  =  14): 


For  — 
r* 


■  41.2  and  0  *  30°,  from  Figure  3(a): 


M'  -  12.75 

(c)  Find  M/M'  from  Figure  3(b)  (Mach  number  trar.sformations-curves ,  7  -  1.4) 

For  M  ■  3.0  and  Q  ■  30° 

e 

M/M'  =  0.922 

(d)  Then: 

M  «  0.922  M'  =  0.922  (12.75) 

M  -  11.75 
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THEORETICAL  NOZZLE  EXHAUST  IN  A  VACUUM  REF  2,  7  =  1.4,  =  2.94 


FIGURE  1  .  THEORETICAL  STREAMLINES  FOR  A  NOZZLE  EXHAUSTING  INTO  A  VACUUM  COMPARED  WITH 
STREAMLINES  FOR  RADIAL  FLOW  FROM  A  POINT  SOURCE 


FIGURE  2.  JET  GEOMETRY  AND  COORDINATE  SYSTEM 


100 


FIGURE  3.  CONTOURS  OF  CONSTANT  MACH  NUMBER  FOR  A  ROCKET  EXHAUSTING 
INTO  A  VACUUM  (BY  PRESENT  APPROXIMATE  SOLUTION) 


M/M’ 


9  (DEGREES) 


R06815 


PIGURE  4b.  MACH  NUMBER  TRANSFORMATION  CURVES 
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FIGURE  5.  A  COMPARISON  OF  THE  APPROXIMATE  AND  CHARACTERISTICS  SOLUTIONS 
FOR  A  ROCKET  EXHAUSTING  INTO  A  VACUUM  (M.  -  2.94,  -y  -  1.4) 


AXIAL  DISTANCE  PARAMETER  (X/re) 


R068  1  6 


FIGURE  6.  COMPARISON  OF  AXIAL  MACH  NUMBER  DISTRIBUTION  FOR  A 
ROCKET  EXHAUSTING  INTO  A  VACUUM  (M..  -  1.0,  1  .4) 


3.93 


O 


FIGURE  7.  A  COMPARISON  OF  THE  APPROXIMATE  AND  CHARACTERISTICS  SOLUTIONS 
FOR  A  ROCKET  EXHAUSTING  INTO  A  VACUUM  (M~  -  3.93,  7  -  1.26) 


FIGURE  8  .  A  COMPARISON  OF  THE  APPROXIMATE  AND  CHARACTERISTICS  SOLUTION  FOR  A  ROCKET 
EXHAUSTING  INTO  A  VACUUM  (M-  -  4 .0 ,  -y  -  1 .24) 


